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Abstract 

We study the dynamical flavour dependence of the lattice esti- 
mates of the heavy meson decay constant in the static limit. We 
perform the analysis by considering negative flavour numbers and by 
extrapolating the results to positive values. We observe a sizeable 
flavour dependence which increases the quenched estimates. 
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The knowledge of leptonic decay constants of heavy mesons enters the 
determination of the Cabibbo-Kobayashi-Maskawa mixing angles. In partic- 
ular, the value of the B meson decay constant constraints the shape of the 
unitarity triangle and sets the size of possible CP violation in B decays. Non 
perturbative estimates of this quantity can be obtained from lattice calcu- 
lations which are still confronted with various sources of systematic errors. 
The presence in the dynamics of a heavy meson system of two very different 
mass scales implies that reasonable values for exploring the B meson dy- 
namics with enough resolution for the heavy quark propagation and enough 
volume for the heavy meson wave function are a lattice size of the order of 
1.5 fm and a lattice spacing of order of 0.02 fm. These values lead to about 
75 lattice points, beyond present computer capabilities. 

Estimates for the physical region of the heavy quark mass are obtained 
by interpolating between the results which are obtained in the charm quark 
mass region and those deriving from an expansion of the fermion action in 
the inverse of the heavy quark mass. The first term in the expansion is the 
so called static limit jl] and corresponds to the approximation where the 
heavy quark does not propagate in space. An important source of possible 
systematic errors of such computations is the quenched approximation: the 
inclusion of dynamical quark loops could alter the quenched estimates and 
vanify the efforts of reducing the statistical fluctuations within the quenched 
approximation. In this letter we make an estimate of the unquenching effects 
on the pseudoscalar decay constant in the static limit with the "bermion" 
method [Q], i.e. we extrapolate from negative to positive flavour numbers the 
results obtained at fixed renormalized parameters. 

On the lattice and in the static limit, the correlator of two local heavy- 
light currents is given by: 

G(t) = ^^yt v , tx+ t(nP^,y)^Sl(x,y)}) (1) 

where S q is the light quark propagator, P(x,y) is the Polyakov line from x 
to y, and L and T are the space and time sizes respectively. 

At large times there is a single state dominating the correlation and 
one can extract the matrix element of the heavy-light current between such 
a state and the vacuum: 
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(2) 



The decay constant in the static limit can be derived from Z^\ 

f P ^M^=V2Z Ren Z L a- 3/2 (3) 

where Z Ren is the suitable renormalization constant of the lattice current. 

In order to simulate theories with negative flavours we introduce a 
"bermion" action which is given by S g [U] + Y,i SfaW, </>%], where S g is the 
standard Wilson action for the gauge sector and the second term is a sum 
over flavour indices of 

s*[^] = £IKM(*)l 2 (4) 

X 

with 4>i(x) the bermion field of flavour i. For the lattice Dirac operator we 
follow the standard Wilson formulation: 

1 1 3 

[QHi x ) = ^-75<Ma:) - 7^75 £ U »i x )( l - ln)M x + ^) 

"^75 £ Ul(x - m)(1 + 7m)&(* - M) (5) 

where k is the Wilson hopping parameter related to the bare mass. 

The square of the Dirac operator is needed for the Montecarlo update 
of the bermion fields and implies that the two point function of these fields 
represents the inverse of such a square. The usual inverse Dirac operator can 
be obtained by remultiplying the bermion field correlations with the Dirac 
operator. 

The correlation G(t) can be estimated with bermion field averages, in 
a way analogous to the one followed in || for the quenched case. In the 
quenched case, the "bermion" fields are not dynamical, they are thermalized 
in a frozen gauge configuration and provide a Monte Carlo inversion of the 
Dirac operator. 

According to the method presented in 0, we define two correlators, 
the first corresponding to eq. ^ 
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1 n b 1 1 

G(t) = -^T,^yt y ,t, + t([Q0^hsP(x,y)^Uy)) (6) 

U b^ J- i= l x ,y ^ 

and a second one corresponding to propagators which are the inverse of the 
square of the Dirac operator: 



1 n b 1 _l_ 

G {Q2 \t) = ^;EE*lAA+»<0l(x)P(x,»)^p^(y)> (7) 

where rif, is the number of bermion flavours over which the sum is performed. 

The statistical fluctuations of these operators are reduced by the "one 
link integral" method M: at fixed bermion configuration one replaces the 
link by its average in the surrounding bermion and gauge configuration. 

As already noticed earlier || the correlation function constructed from 
the inverse of the square of the Dirac operator in eq. |7| projects more preco- 
ciously on the lowest lying state. We have therefore extracted the energy shift 
AE from this correlation and then fixed its value in the fit of the canonical 
correlator of eq. |6|. The simulations were performed on a 25 Gigaflop machine 
of the APE series. The update procedure was for the gauge sector a Cabibbo- 
Marinari pseudo-heatbath || followed by three overrelaxation sweeps and 
for the bermions a heat bath followed by ten overrelaxation sweeps . The 
measurements were performed every five sweeps: we have checked that the 
autocorrelation time, defined as r = crf rue /(2er^ ai y), where a true is evaluated 
by collecting measurements into bins and o n dif is evaluated by considering 
each measurement as an independent one, is around 0.5 for the correlations 
G(t) and around 1 for \t). According to the procedure of extrapolating 
at fixed renormalized quantities, the values of the (3 and k in each simulation 
with a given number of bermions were chosen to match the p and 7r masses 
at the corresponding quenched values. We have explored two values of the 
ratio i?2 = ^/ m p) 0-5 an d 0.6, with a lattice spacing corresponding to the 
value (3 = 5.7 of the quenched case and appropriate values of k. 

We have made simulations with nf = —2, —4 and —6 for R 2 = 0.5 and 
rif = —2, —4, —6 and —8 for R 2 = 0.6. The quenched values are those of ref. 
H with a slight different statistics in the case of k = 0.163. The errors are 
evaluated as follows: we form clusters of 200 measurements. On these clusters 
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we apply a jackknife algorithm to estimate the errors of the correlations. Each 
"jackknife cluster", defined by single elimination of a cluster from the total 
sample is fitted with the MINUIT minimization program. From the spread 
of the fits on each jackknife cluster we extract the errors that we quote for 
our results. For the operator of eq. we have checked that the results of 
one and two mass fits are generally compatible while for the operator of eq. 
H a two mass fit is necessary with the lowest mass fixed by the fit of the 
other operator. We quote separately the statistical error of the fit and the 
error deriving from the uncertainty by which the correlation \t) fixes 
the value of the lowest energy shift. All the results have been obtained on a 
16 3 x 32 lattice and the correlations have been taken in all directions, with a 
caveat. The correlations in the space directions, for values of the separation 
close to the boundary (beyond 11-12) may suffer from finite size effects due to 
the presence of dynamical "bermions" . Indeed, the light quark can connect a 
point in the original lattice with a point sitting in the backward copy, while 
the heavy quark is explicitely propagated only in the forward direction and 
sits always in the original lattice. The result is a path wrapping around 
the lattice. In the quenched approximation and in the confined phase, these 
paths average to zero, but they are in general non zero in the unquenched 
case. The fermion determinant may create a fermion loop wrapping around 
the lattice which, through its correlation with the one from the heavy-light 
propagator can give to the latter a non zero expectation value. The effect 
is visible beyond the statistical errors for t > 10. We have then included in 
the fits the points from the correlation in the space directions up to t = 10 
only. For the time direction which is twice longer the effect is negligible for 
the maximum values (t = 14) of the time distance where the signal is still 
appreciable. 

Figures 1 and 2 give the flavour dependence of the static decay con- 
stants for the two values of B,2- in Tables 1 and 2 we present our detailed 
numerical results for the energy shift, for and for the parameters of our 
simulations including the matched values of rr and p masses. The extrap- 
olation of the results to n/ = +3 has been attempted with a linear and a 
quadratic fit with the results collected in table 3. While the former varies 
with the number of negative flavours included in the fit, the latter is more 
stable against this variation and leads to values which are compatible within 
our errors. This applies also to the case of R2 = 0.6 where we have explored 
the flavour dependence up to nj = —8. In order to take into account the 
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error arising from an imperfect matching of some simulations, we have ex- 
trapolated the values of Z^ corrected by the factor (m^ - /m^) 3 / 2 . The 
unquenched estimates at R 2 = 0.5 and R 2 = 0.6 can be further extrapolated 
in the quark mass (i.e. in the inverse of the Wilson hopping parameter) to 
the chiral limit where a physical values of f^ attc can be obtained by using 
the experimental value of the p mass as an input. 

Our final results for the unquenched f^ attc at a lattice spacing corre- 
sponding to (3 = 5.7 of the quenched case, extrapolated to the chiral limit 
with n/ = 3 and not corrected for the renormalization constant is: 

/ o m P h y s 

f ^ c/z n en = Zl] J_ (_4_)3/ 2 = .66(7)GeV 
to be compared with the quenched result: 

f static I z Ren = ().50(3)GeV 

The unquenched estimate for the ratio f§^ tlc / f 'B^ Uc cannot be obtained 
directly from the quark mass dependence of our results. Indeed, in our sim- 
ulations we have identified the valence and the sea quark masses. This is 
appropriate for the light quark spectrum, while for the strange quark spec- 
troscopy one should keep the two types of masses distinct. We defer such a 
refined study to a future work. 

Including dynamical fermion loops increases the static decay constants, 
an effect already claimed by the MILC collaboration which uses dynamical 
staggered quarks 0. 

The results discussed so far are for the bare quantities, not corrected 
for the appropriate renormalization constants. Their inclusion partially com- 
pensates the observed flavour dependence, but the effect can be estimated 
to be smaller than the increase of Zl due to the unquenching. Indeed by 
using the bare couplings in the one loop expression for Z Ren , the percentage 
of the unquenching effects is decreased by a few units. By using an improved 
perturbation theory, i.e. a coupling constant in a Parisi || scheme, the effect 
is larger and the percentage decreases by about ten units. The extrapolated 
values of Z Ren in the Parisi coupling scheme are given in Table 3. 

A sizeable flavour dependence survives the corrections of the renormal- 
ization constant, quantifies into a 10-15 % increase in the intermediate quark 
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mass region that we have explored and becomes about 20% in the chiral 
limit. 
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n f 




P 


K 


AE 


Z L 




m p 


-8 


6000 


6.99 


0.153 


0.742(2) 


0.440(2) [7] 


0.563(5) 


0.712(10) 


-6 


5000 


6.71 


0.1548 


0.751(3) 


0.466(2) [8] 


0.559(6) 


0.716(9) 


-4 


3000 


6.4 


0.157 


0.770(3) 


0.515(4) [9] 


0.552(3) 


0.710(6) 


-2 


4400 


6.06 


0.1598 


0.789(2) 


0.573(2) [7] 


0.566(4) 


0.714(8) 





9200 


5.7 


0.163 


0.811(1) 


0.648(4) [4] 


0.562(2) 


0.719(4) 



Table 1: Results for AE and Zl at R2 = 0.6 and the matched values of 
m w and m p . The second column gives the number of measurements in each 
case. For Zl we separate the statistical error () from the one coming from 
the determination of AE []. 



n f 




P 


K 


AE 


z L 


m w 


m p 


-6 


3600 


6.785 


0.156 


0.733(2) 


0.440(2) [5] 


0.477(8) 


0.67(1) 


-4 


3400 


6.463 


0.158 


0.744(2) 


0.465(7) [6] 


0.464(4) 


0.658(8) 


-2 


5000 


6.1 


0.161 


0.762(3) 


0.520(10) [9] 


0.467(2) 


0.662(5) 





6200 


5.7 


0.165 


0.786(2) 


0.590(7) [7] 


0.457(3) 


0.659(8) 



Table 2: The same quantities of table 1 for R 2 = 0.5. 
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#of 
points 


#2 = 0.5 


Z L (n f = 3) 


Z Ren (n f = 3) 


Z L Z Ren (n f = 3) 


Z L Z uen (n f =3) 


Z L Z Ren {n f =0) 


parabolic 


linear 


parabolic 


linear 


parabolic 


linear 


parabolic 


2 
3 
4 


0.744(55) 
0.730(23) 


0.702(38) 
0.681(22) 
0.655(17) 


0.647 
0.649 


0.662 
0.670 
0.678 


0.490(42) 
0.485(18) 


0.473(29) 
0.465(17) 
0.453(14) 


1.16(4) 




K2 = l 


.6 


Z L (n f = 3) 


Z Ren (n f = 3) 


Z L Z Ren (n f = 3) 


Z L Z Hen (n f =3) 
Z L Z Be "(n f =0) 


parabolic 


linear 


parabolic 


linear 


parabolic 


linear 


parabolic 


2 
3 
4 
5 


0.779(42) 
0.760(19) 
0.778(12) 


0.742(25) 
0.730(17) 
0.730(13) 
0.718(12) 


0.659 
0.659 
0.661 


0.668 
0.673 
0.678 
0.684 


0.519(33) 
0.506(15) 
0.522(9) 


0.503(20) 
0.499(13) 
0.504(11) 
0.498(9) 


1.14(2) 



Table 3: We give the results for the values of Zl extrapolated to rif = 3 at 
R 2 = 0.5 and R 2 = 0.6 as a function of the number of points used in the 
extrapolation. 
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n f 

Figure 1: The values of Z L as a function of the flavour number for R 2 = 0.5. 
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n f 

Figure 2: The values of Z L as a function of the flavour number for R 2 = 0.6. 
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